Quasi One-Dimensional Bosons in Three-dimensional Traps: 
From Strong Coupling to Weak Coupling Regime 



L. Salasnich 1 , A. Parola 2 and L. Reatto 1 
1 Dipartimento di Fisica and INFM, Universitd di Milano, 
Via Celoria 16, 20133 Milano, Italy 
2 Dipartimento di Scienze Fisiche and INFM, Universitd dell'Insubria, 
Via Valeggio 11, 23100 Corno, Italy 

We analyze a recent experiment on a Tonks-Girardeau gas of 87 Rb atoms (T. Kinoshita, T. Wenger, 
and D.S. Weiss, Science 305, 1125 (2004)). We find that the experimental data are compatible with 
the one-dimensional theory of Lieb, Seiringer and Yngvason (Phys. Rev. Lett. 91, 150401 (2003)) 
but are better described by a theory that takes into account variations in the transverse width of 
the atomic cloud. By using this theory we investigate also the free axial expansion of the 87 Rb 
gas in different regimes: Tonks-Girardeau gas, one-dimensional Bose- Einstein condensate and three- 
dimensional Bose-Einstein condensate. 
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Two experimental groups [1,2] have reported the ob- 
servation of the one-dimensional (ID) Tonks-Girardeau 
(TG) gas [3,4] with ultracold 87 Rb atoms in highly 
elongated traps. A rigorous theoretical analysis of the 
ground-state properties of a uniform ID Bose gas, in- 
cluding the beyond-mean-field TG regime of impenetra- 
ble bosons, was performed by Lieb and Liniger (LL) forty 
years ago [4]. Recently, motivated by the experimental 
achievements, an extension of the LL theory for finite and 
inhomogeneous ID Bose gases under longitudinal con- 
finement has been proposed on the basis the local density 
approximation (LDA) [5]. Ohberg and Santos [6] have 
suggested that the LDA is improved by including a gra- 
dient term that represent additional kinetic energy asso- 
ciated with the inhomogeneity of the gas. This conjecture 
has been rigorously proved by Lieb, Seiringer and Yng- 
vason [7]. More recently we have introduced [8] a vari- 
ational approach, called generalized Lieb-Liniger theory 
(GLLT), which reduces to the Lieb-Seiringer- Yngvason 
theory (LSYT) in the ID regime and, in addition, gives 
an accurate description of the crossover from the ID 
regime to 3D regime. 

In this Brief Report we apply the GLLT to analyze 
the experimental results of Kinoshita, Wenger and Weiss 
[1] . Contrary to the data of Paredes et al. [2] which are 
deeply in the TG regime, the data of Rcf. [1] cover dif- 
ferent quantum-dimensional regimes and we show that 
the experimental atomic cloud is better described by the 
GLLT than the LSYT. The GLLT is then used to de- 
termine the axial free expansion of the 87 Rb gas. We 
predict a self-similar expansion whose growth is strictly 
related to the quantum-dimensional regime of the initial 
configuration. 

In the LSYT [7] the longitudinal density p(z) of a zero- 
temperature Bose gas is obtained by minimizing the fol- 
lowing energy functional 
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where g = 2a s /a 2 L is the interaction parameter, with a s 
the s-wave scattering length and aj_ = (h/(mu>±) 1 ' 2 the 
characteristic length of the transverse harmonic poten- 
tial with frequency uj±, V(z) is the longitudinal external 
potential and e(x) is the Lieb-Liniger function, which 
is defined as the solution of a Fredholm equation and 
it is such that e(x) w x — 4/(37r)x 3 / 2 for x <C 1 and 
e(x) « {tt 2 /3)(x/(x + 2)) 2 for x > 1 [4]. 

As previously stressed, the LSYT is valid in the pure 
ID regime, where the transverse width R± of the Bose 
gas is frozen and equal to the harmonic length a±_ . In our 
GLLT [8] the transverse properties of the Bose gas are 
taken into account by considering the adimensional width 
a = R±/a± as a variational parameter, i.e. a = <t(z). 
The theory gives the following energy functional 
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where Huj±(<t~ 2 + a 2 )p/2 is the transverse energy den- 
sity of the Bose gas. Note that in the low density re- 
gion where a ps 1 (ID regime) the functional -Egll re- 
duces to Elsy- Instead, at higher densities (3D regime) 
where e(g/(pa 2 )) « g/(pa 2 ), the functional Eqll gives 
the energy functional of the nonpolynomial Schrodinger 
equation (NPSE) [9], an effective ID differential equation 
we have derived from the 3D Gross-Pitaevskii equation 
[10] to describe Bose-Einstein condensates (BECs) under 
transverse harmonic confinement. 

Taking into account the normalization condition of the 
longitudinal density p(z) to the total number iV of atoms, 
the minimization of the energy functional Eqll with re- 
spect to p(z) gives the equation 
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where p, is the chemical potential fixed by the normal- 
ization condition. The minimization of the energy func- 
tional Eqll with respect to a(z) gives instead the equa- 
tion 
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This implicit equation must be solved numerically but 
analytical results can be obtained in limiting cases. In 
particular, under the condition a s <C a± [11], one gets 
a ps y/al(pg) 1 ' A for p > l/a s (3D BEC regime) and 
er ps 1 for p <C X/a s (ID regime). Note that the ID 
regime contains two subregimes: the ID BEC regime 
for a s ja\ <C p <C l/a s , and the TG regime, where 
e(g/(pa 2 )) ps tt 2 /3, for p < a s /a^_. 
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FIG. 1 (color online). Axial density p(z) of a Bose gas of 
50 87 Rb atoms in an anisotropic harmonic trap with frequen- 
cies oj± and uj z . Four values of trap anisotropy: r = \/\ ma x, 
where A = lu±/lu z and \ m ax = 2570. Comparison among 
various theories: GLLT (full line); GLLT with TF approx- 
imation (dotted line); LSYT (dashed line); LSYT with TF 
approximation (long dashed line). Length in units a z , and 
axial density in units l/a z . 

In the experiment of Kinoshita, Wenger and Weiss [1] 
an ensemble of about 6400 parallel ID traps has been 
created by means of a 2D optical lattice, that strongly 
confines atoms in ID tubes. The traps differ only in the 
number of 87 Rb atoms each contains. As discussed in [1], 
each trap can be modelled by an anisotropic harmonic 
potential, where the axial frequency is co z = 2it x 27.5 Hz 
and the transverse frequency u>j_ is tuned by changing 
the energy depth f/ of the confining optical lattice. The 
maximum transverse frequency is 2ir x 70.7 kHz; it fol- 
lows that the maximum trap anisotropy is X m ax — 2570, 
where A = lu±/lo z . Kinoshita et al. [1] have worked with 
r = \/\„iax ranging between r — 1 and r = 10 _1 , a 
regime where the Bose gas is quasi-lD. 



By adding an imaginary time term at right side of 
Eq. (3) and solving Eq. (3) and Eq. (4) with a finite- 
difference Crank-Nicolson predictor-corrector scheme 
[12], we calculate the density profile p(z) of the Bose 
gas using the trap parameters of Ref. [1]. In Eq. (3) we 
set V(z) = m,L0 2 z z 2 /2. The results are shown in Fig. 1, 
where we compare LSYT with GLLT, plotting also their 
Thomas-Fermi (TF) approximations which neglect the 
gradient term in Eq. (1) and Eq. (2). Figure 1 shows 
that, for the trap anisotropy r = 1, the TF approxima- 
tion is remarkably accurate but there are instead quan- 
titative differences between LSYT and GLLT due to the 
fact that a(p(z)) is not strictly equal to one. The differ- 
ences between LSYT and GLLT increase by reducing the 
anisotropy r of the trap, and also the TF approximation 
worsens. In fact, by reducing r one induces a crossover 
from the ID regime to the 3D regime, whose description 
cannot be accounted for by LSYT that is purely ID [7]. 
We have verified that the differences between LSYT and 
GLLT are instead reduced by taking larger values of the 
anisotropy ratio r. 
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FIG. 2 (color online). Root mean square (RMS) full 
length of the Bose cloud of 87 Rb atoms. Uo/E rec is the 
trap depth Uo in units of the atom recoil energy E rec , where 
Uo/Erec = 87(A/A mol ) 2 with A = lu±/lu z and \ max = 2570. 
Experimental data are obtained in [1] by averaging over sev- 
eral clouds with different numbers of atoms (maximum num- 
ber Nmax = 54) . 

In Fig. 2 we plot the experimental data of the root 
mean square (RMS) length of the atomic cloud obtained 
by Kinoshita, Wenger and Weiss [1]. This RMS axial 
length has been derived as the average over several clouds 
which have different numbers of atoms. In Fig. 2 the 
RMS full length is plotted versus U /E rec = 87r 2 , where 
Uq is the energy depth of the optical potential and E rec 
is the recoil energy of the gas. Kinoshita et al. [1] have 
compared their data with LSYT [7] in the TF approxima- 
tion [5] . By performing the averaging procedure of Ref. 
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[1] on our LSYT results, we obtain the same theoretical 
curve shown in Fig. 4 of Ref. [1] and plot it (dashed 
line) in our Fig. 2. With the same averaging procedure 
we obtain also the GLLT curve (dotted line). Figure 2 
shows that our GLLT fits better than LSYT the experi- 
mental data. As discussed in [1], the experimental data 
for Uo/E rec < 20 are strongly affected by the tunnel- 
ing between adjacent tubes and the theoretical analysis 
based on a single tube overestimates the axial width of 
the Bose cloud. 

Our GLLT is a simple and useful tool to investigate 
also dynamical properties of a Bose gas under transverse 
harmonic confinement. The dynamics of the Bose gas can 
be described by means of a complex classical field &(z, t) 
that satisfies the equation ifi <9t$ = j^Egll[\&\ 2 , 
where \<&(z, t)\ 2 — p(z, t). This time-dependent nonlinear 
Schrodingcr equation can be explicitly written as 
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where p[p] is the bulk chemical potential, given by 
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Here we analyze the free axial expansion of the Bose 
gas by means of the time-dependent GLLT. We con- 
sider the expansion of the Bose gas when the ax- 
ial harmonic confinement is removed (lu z = 0) while 
the radial one is kept fixed. By setting $>(z,t) — 
\J p(z, t) exp (iS(z, t)/K), the Eq. (5) is equivalent to 
the two hydrodynamics equations d t p + d z (pv) = and 
m d t v + d z [-Ti 2 /(2m^/p)d 2 ^p + mv 2 /2 + V + p] = of a 
ID viscousless fluid with density field p(z, t) and velocity 
field v(z, t) = (H/m)d z S(z, t). If the initial axial width of 
the cloud is larger than the healing length h/\/2mp then 
one can safely neglect the quantum pressure (QP) term 
—h 2 /(2my r p)d zy /p in the second hydrodynamics equa- 
tion. Note that the bulk chemical potential p of Eq. (6) 
scales as p 1 in the three relevant regimes: 7 = 2 in the 
TG regime, 7 = 1 in the ID BEC regime, and 7 = 1/2 in 
the 3D BEC regime. It is straightforward then to prove 
that in these regimes the cloud density decreases dur- 
ing the time evolution following the self-similar solution 
p(z,t) — p(z/b(t),t = 0)/b(t), where the adimensional 
axial width b(t) satisfies the equation 
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The Eq. (5) with Eq. (6) must be solved self-consistently 
with Eq. (4). It is important to observe that for p <C 
°s/ a i (TG regime) one has a s=s 1 and p s=s ft 2 it 2 p 2 / (2m) 
and Eq. (5) reduces the the time-dependent nonlinear 
Schrodinger equation (NLSE) introduced by Kolomieski 
et al. [13]. For p <C l/a s (ID regime or TG regime) 
still we have a w 1 but p w 3ft 2 p 2 e(g / p) / (2m) — 
h 2 gpe'(g/p)/(2m) and Eq. (5) becomes equal to the 
time-dependent NLSE proposed by Ohberg and Santos 
[6], which reduces to that of Kolomieski et al. in the 
TG regime. For p 3> a s /a^ (ID regime or 3D regime) 
the Eq. (5) gives a w (1 + ga^p) 1 / 4 , the Eq. (6) reads 
p w h 2 gp/(ma 2 ) + Huj ± (a' 2 + a 2 - 2)/2 and Eq. (5) 
becomes the time-dependent NPSE [9]. The NPSE coin- 
cides with the equation of Ohberg and Santos in the ID 
BEC regime, where a s ja\ <C p <SC l/a s , but it accurately 
describes also the 3D BEC regime, where /)> l/a s and 

a w V^I(5P) 1/4 - 

The time-dependent (TD) GLLT given by Eq. (5) 
must be carefully employed. In fact, Girardeau and 
Wright [14] have shown, studying the interference of 
two Bosonic clouds, that the time-dependent NPSE of 
Kolomieski et al. [13] overestimates the coherence of the 
Tonks-Girardeau gas. However, the TD GLLT can be 
safely used to calculate collective properties of the Bose 
gas. In Ref. [8] we have shown that the time-dependent 
GLLT gives the expected values for the axial breathing 
mode J7 Z of the Bose cloud in a harmonic confinement 
with V(z) = muj 2 z 2 /2. In particular, we have found 
VL Z = 2Q Z in the TG regime, Sl z = y/3uj z in the ID BEC 
regime, and Q z = y/5/2uj z in the 3D BEC regime. 



The solution of this equation with initial conditions 
6(0) = 1 and 6(0) = can be obtained by quadratures. 
For large t the solution is b(t) = \j2j^ u> z t. In general, 
p(p) is not a power law and during the expansion the 
functional dependence of p(p) changes throughout the 
whole cloud: the expansion is no more truly self-similar. 
However, by evaluating the dynamics at the center of the 
cloud (z = 0) where the initial density is po, from the two 
hydrodynamics equations without the QP term one finds 
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For large t the solution of this equation is b(t) = 
\/2/7 Lu z t, where naturally appears the effective poly- 
tropic index 
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that is the logarithmic derivative of the chemical poten- 
tial p. In order to verify these analytical results we sim- 
ulate the free axial expansion of the Bose gas by numer- 
ically solving Eq. (5) self-consistently with Eqs. (4,6). 
In our simulation we employ a real-time Crank-Nicholson 
method [12]. 

In Fig. 3 we plot the time evolution of the normalized 
axial width W(t) of the Bose gas with N = 50 atoms 
and initial anisotropy of the harmonic trap given by r = 1 
(luj_/lu z = 2570). After a transient the width W(t) grows 
linearly in time, as clearly shown by the lower panel of 
Fig. 3, where we plot the width velocity dW/dt. In 
Fig. 3 we plot also W(t) and dW/dt calculated with two 
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alternative methods: by using Eq. (8) and by using Eq. 
(7) with 7 given by Eq. (9), in both cases taking into 
account Eqs. (4,6). 




FIG. 3 (color online). Axial free expansion of the Bose 
gas of 50 87 Rb atoms with initial anisotropy lu±/lo z = 2570. 
Upper pa nel: time evolut ion of the normalized axial width 
W(t) — yl (z 2 }(t) I (z 2 )(0) of the atomic cloud. Lower panel: 
time evolution of the velocity dW/dt of the axial width. TD 
GLLT given by Eq. (5) and numerical solution of Eq. (7) 
with 7 from Eq. (9). 




FIG. 4 (color online). Slope a = y2/7 of the asymptotic 
linear growth of the axial width of the Bose gas as a function 
of the initial density po of the 87 Rb cloud at the center of the 
trap. Three values of the initial anisotropy of the harmonic 
trap: anisotropy ratio r defined as in Fig. 1. 

L.S. thanks D.S. Weiss for enlightening e-suggestions. 



Figure 3 shows that the approximation based on the 
polytropic index is reliable; note that in all our tests the 
relative error in the slope of the asymptotic linear growth 
is always within 5%. 

In Fig. 4 we plot the slope a = y/2/'y as a function 
of the initial density po at the center of the cloud, a 
smoothly changes from 1 to \J~2~ and than to 2 in the 
transition from the TG regime to the ID BEC regime 
and then to the 3D BEC regime. Figure 4 shows that for 
r <C 1 in addition to the two plateaus corresponding to 
a = 2 at high densities (3D BEC regime) and to a = 1 at 
low densities (TG regime), an additional plateau appears 
with a — y2 at intermediate densities (ID BEC regime). 

In conclusion, we have analyzed a recent experiment on 
dilute Bosons in the Tonks- Girardeau regime by using a 
beyond-mean-held theory that takes into account trans- 
verse variations of the atomic cloud. Remarkably, our 
theory shows a better agreement with the experimental 
results than the simple ID approach proposed by various 
authors in the past years. On the basis of this theory, 
we have predicted a quasi self-similar expansion of the 
gas when the axial confinement is removed, showing that 
asymptotic linear growth of the axial width strongly de- 
pends on the quantum-dimensional regime of the Bosonic 
cloud. 
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